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Abstract 

The discrete non-linear Schrodinger (NLS) model in the presence of an integrable defect 
is examined. The problem is viewed from a purely algebraic point of view, starting from 
the fundamental algebraic relations that rule the model. The first charges in involution 
are explicitly constructed, as well as the corresponding Lax pairs. These lead to sets of 
difference equations, which include particular terms corresponding to the impurity point. A 
first glimpse regarding the corresponding continuum limit is also provided. 
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1 Introduction 

The presence of defects in 1 + 1 integrable field theories has been the subject of intense 
research during the recent years (see e.g. [1] [11]). It is well established by now that 
the requirement of integrality leads to a set of severe algebraic constraints that should be 
satisfied by the associated degrees of freedom as well as the relevant physical quantities, 
such as scattering matrices, at the quantum level (see e.g. [HI2])- In integrable field theories 
the defect is usually introduced as a discontinuity together with suitable sewing conditions 
[3] [9]. In this case no systematic algebraic description exist so far with the exception of 
some recent attempts [9], but again the issue of integrability is not fully resolved to our 
understanding (see also [1 Oj ) . 

In the present study we start our investigation using an integrable model on the one 
dimensional lattice, that is the discrete non-linear Schrodinger (NLS) model, and impose an 
ultra-local integrable defect (see also relevant considerations in |11]). We by construction 
deal with an integrable system by a priori imposing the necessary algebraic constraints that 
ensure the integrability of the model. We then explicitly construct the first integrals of 
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motion as well as the relevant Lax pairs. The corresponding equations of motion are also 
derived. The pertinent question in this frame is whether and how integrability is preserved 
in the continuum limit. In other words is the underlying algebra that ensures integrability 
modified and how? We make a first attempt to answer these questions by considering the 
continuum limits of the first couple of integrals of motion as well as the corresponding Lax 
pairs; then certain continuity or sewing conditions are naturally induced via the process. 
These preliminarily results on the continuum limit, provide a first insight on how a systematic 
continuum process should be formulated within this context. 

The outline of the present article is as follows: In the next section we briefly review 
the discrete NLS model with periodic boundary conditions. In section 3 the first three 
integrals of motion are derived together with the corresponding Lax pairs. We also derive 
the equations of motion associated to the third charge, which is the typical Hamiltonian. 
In the next section we consider the discrete NLS model in the presence of an integrable 
local defect. We introduce a suitable Lax operator associated to the defect point so that 
the system is by construction integrable. We then derive the first integrals of motion as 
well as the corresponding Lax pairs. These are novel expressions that contain non-trivial 
contributions due to the presence of the integrable defect. The related equations of motion 
for the Hamiltonian are also derived. Particular emphasis is given exactly on the defect 
point where the equations of motion are of a completely different form compared to the 
other points due to the structural dissimilarity of the associated Lax operator. Finally, in 
section 5 we provide a first glimpse on the continuum limits of the derived physical quantities. 
Certain sewing conditions naturally arise as continuity requirements, ensuring the Poisson 
commutativity of the first two continuum charges. 



2 The periodic DNLS model 

We shall briefly review the discrete NLS model with periodic boundary conditions. We 
shall reproduce the first three local integrals of motion and the associated Lax pairs. In 
the subsequent sections we shall repeat these derivations in the presence of an integrable 
defect. Our ultimate aim is to make a contact with recent results on defects arising in 1+1 
integrable field theories (see e.g. [5]). We aim at taking the appropriate continuum limit, 
that will provide the classical continuum models with defects that are still integrable. This 
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is conceptually and technically a very intriguing problem, and will be pursued in full detail 
in forthcoming investigations. Here, however we provide some preliminary results regarding 
the continuum limit of the first two integrals of motion. 



3 Local Integrals of motion 

Our main aim in this section is to extract the first integrals of motion for the periodic discrete 
NLS model. The associated Lax operator is given by (see e.g. [12]): 



L aj (\) = XD.+A, 

nr. ■ \ 

(3.1) 



X + Nj xj 
-X, 1 



where Nj = 1 — XjXj. 

We shall focus here, mainly for simplicity, on the classical case. However, we have to 
note that our results are valid in the quantum well. The L matrix satisfies the 



fundamental algebraic relation 



{i a (Ai), L h {\ 2 )} = [ ra6 (Ai - A 2 ), L a (Ai) L h {X 2 )\ . (3.3) 
In this case the r-matrix is the familiar sl 2 Yangian matrix [14] : 

r(A) = jV, (3.4) 

V is the permutation operator: v(E®b^ = 6® a. The formula (13. 3 j) is then realized by the 
following relations: 

{xi, Xj} = Sij, (3.5) 

that is x, X are canonical variables. 

The discrete model with N sites, and periodic boundary conditions is associated to the 
transfer matrix defined as [T3], [151 EH] : 

t{\) = Tr a T a {\) where T a (X) = L aN {X)L aN ^(X) . . . L al (X). (3.6) 



1 The same L-operator holds for the quantum case as well. It then satisfies: 

JZia(Ai - A a ) £i(Ax) L 2 (X 2 ) = L 2 (A 2 ) L^Ai) R 12 (\i - A 2 ), (3.2) 
where R(X) = AI + V. 



T is the monodromy matrix also satisfying the quadratic algebraic relation (13. 3p . In the 
notation L 0i , the index a denotes the auxiliary space, whereas the index i denotes the ith 
site on the one dimensional lattice derived by (13. 6p . As will be transparent later in the text 
in the continuum limit the discrete index i will be replaced by the continuum coordinate x. 

The transfer matrix t(X) as is well known provides all the charges in involution. Indeed, 
via (13. 3 p one readily shows that 

{i(A), t(Ai)} = 0, (3.7) 

hence the system is by construction integrable. 

To derive the local integrals of motion one should expand the lnt(A) in powers of A or ~. 
In this case we expand in powers of j-, because the L-matrix (13. ip reduces to the degenerate 
matrix D at A — > oo. 

Let us now expand the monodromy matrix: 

1 - 

T(A^oo) ex D N ...D 1 + -J2D N ...D i+1 A i D i _ 1 ...D 1 

A .=1 

+ -Ja Yl Dn ■ ■ ■ A+i4A-i • • • D j+1 Aj ...D 1 

i>j 

+ ^ Yl Dn ■ ■ ■ A+iAA-i • • • D j+1 Aj . . . D k+l A k ...D 1 

i>j>k 

+ ... (3.8) 

Taking into account the latter expansion and the definition of the transfer matrix we con- 
clude: 

\nt(X^oo)ocjH 1 + ^H 2 + ^H 3 + ..., (3.9) 
where the extracted integrals of motion have the following familiar form (see also e.g. [i~2|[T7]) 

N 



i=l 

N . N 

# 2 = - ]Twq - - n: 



2 

i=l i=l 



N N N 



3 

i=l i=l i=l 



The latter provide the first integrals of motion (number of particles, momentum and Hamil- 
tonian respectively) of the whole hierarchy for the NLS model. It is clear that the continuum 
limits of the above quantities provide the corresponding integrals of motion of the continuum 
NLS model p21 EJ- The latter expressions are valid in the quantum case as well (see e.g. 

punj). 

3.1 The Lax pair formulation 

Let us now briefly review how the Lax pair associated to each local integral of motion is 
derived via the r- matrix formulation (see also [IS])- Introduce first the Lax pair (L, A) 
for discrete integrable models, and the associated discrete auxiliary linear problem (see e.g. 

ipj = Aj ipj. (3.11) 

From the latter equations one may immediately obtain the discrete zero curvature condition 
as a compatibility condition: 

Lj = Aj + i Lj — Lj Aj. (3.12) 

Recall that the index j denotes the site on an one dimensional lattice, and it will be replaced 
in the continuum limit by the continuum coordinate x. In the continuum limit as will be 
clear the equations (13. lip . ( I3.12p reduce to the continuum linear auxiliary problem and the 
continuum zero curvature condition respectively (see section 5). 

Let us introduce at this point some useful notation. We define for i > j: 

T a {i,j;X) = L ai {X)L ai - 1 {X)...L aj {X). (3.13) 

To be able to construct the Lax pair we should first formulate the following Poisson structure 
H3|: 

{r a (A), L bi Cu)} = T a {N,j + l;X)r ab (X-»)T a {jA;X)L bj {f Jl ) 

- L bj (fi)T a (N,j;X)r ab (X- fj,)T a (j -1,1;X). (3.14) 
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It then immediately follows for the generating function of the local integrals of motion: 

{ lnt(A), L bj -(/x)} = r x (A) Tr a (T a {N,j + 1; A) r ab (\ - T a (j, 1; A)) L fej (/i) 

- Ly(^) r x (A) Tr a (T a (N,j;X) r ab (X-fi) T a (j - 1, 1; A)).(3.15) 

Recalling the classical equation of motion 

L,(/i) = {lnt(A), L,(/i)}, (3.16) 

and comparing with expression (13.151) we obtain 

A J -(A,A*)=r 1 (A) tr a [r a (iV,j;A) r ab (A-/i) T a (j - 1, 1; A)j , (3.17) 

where the relevant classical r-matrix is given in (13. 4p . 

Substituting the r-matrix into the latter expression we conclude that 

A, (A, n) = t -^- T(j - 1, 1; A) T(N,j; A). (3.18) 
A — n 

Expansion of the latter expression in powers of j provides the Lax pairs associated to each 
one of the local integrals of motion (see also [IE]), i-e.: 



0/ J \ -Xj-i 



(3) / (1 + XjXj_i flXj — XjNj + Xj + i 



j 



(3.19) 



Both the Lax pair via the zero curvature condition and the Hamiltonian description 
give rise to the same equations of motion. Consider for instance the equations of motion 
associated to H 3 (and the Lax pair L, A^). Indeed from 

Xj = {H 3 , Xj }, X 3 = {H 3 , X,}, (3.20) 

and via the zero curvature condition for the pair L, we obtain the following set of 
difference equations: 

xj = x j+2 - 2x j+1 Nj - Xj+xNj+i + XjNj + x 2 jXj_i + x j+1 
Xj = -X j ^ 2 + 2X 3 _ 1 N 3 +X^ 1 N^ 1 -X J N 2 J -X]x j+1 -X :j _ 1 . (3.21) 

With this we conclude our brief review on the periodic discrete NLS model. 
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4 The DNLS model with integrable defect 




We shall henceforth focus on the discrete NLS model in the presence of an integrable defect. 
We shall basically extract the local integrals of motion and the corresponding Lax pairs for 
the aforementioned model, and shall derive the modified equations of motion due to the 
presence of the defect. 

Let us first describe the algebraic setting for the defect per se. Introduce the Lax operator 
associated to the defect, which is located at a particular site say n: 

La., 

X ) (1.1) 

the index n simply denotes the position of the defect on the one dimensional spin chain. 
Note that the L matrix is required to obey the same Poisson bracket structure with the bulk 
matrices L ( 13. ip so that integrability is ensured. The entries of the above L matrix may be 
parameterized as (see e.g. [HI], and references therein) 

a n = -6 n = ^cos(20 n ), /? n = isin(2# n )e 2 ^, ln = ± sin(2# n ) e - 2 ^, (4.2) 

the fields 9 n , 4> n may be rewritten in terms of the canonical variables p n and q n as 

C0s(26y = p n , (p n = q n 

tin, Pn} = t- (4.3) 

It is then immediately shown via the algebraic relation (13. 3p that the elements a n , (3 n , j n , 5 n 
satisfy the following exchange relations: 

{«„, 7„} = -7n 

{P n , ln } = 2a n (4.4) 
which are the typical sl2 exchange relations. 

For simplicity, and in order to avoid unwanted boundary effects we shall consider the 
defect away from the ends of the one dimensional lattice model. Inserting the defect at the 
n site of the one dimensional lattice the corresponding monodromy matrix is expressed as: 

T a (A) = L aN (X)L aN ^(X) . . . L an (X) . . . L a i (A). (4.5) 



7 



Note that due to the fact that the L-operator is required to satisfy the same fundamental 
algebraic relation as the monodromy matrix, the trace of it -the transfer matrix- provides 
a family of Poisson commuting operators. Having the latter expression at our disposal we 
may now construct the desired physical quantities. 

4.1 Local Integrals of motion 

First we wish to extract the associated local integrals of motion. They are obtained, as in 
the previous section, from the expansion of lnt(A). Let us first present the expansion of the 
relevant monodromy matrix: 

T(X ->■ oo) oc D N ... D 1 
+ 



\(Y1 Dn --- A+i^i A-i ■■■D 1 + D N ... D n+x A n D n _ x ...DA 
+ ^J2 Dn --- A+i^A-i • • • D j+1 Aj ...D 1 

i>j 

+ ^ Yl Dn ■ ■ ■ Dn+lKDn-! . . . D j+1 Aj . . . D 1 

n>j 

+ ^ DN ■ ■ ■ D j+1 A 3 D 3-1 ■ ■ ■ D n+lAn ■ ■ ■ D X 

j>n 

+ ... (4.6) 
The technical details are omitted for brevity, and we directly provide the final expressions: 

log t(X) = iftx + ±H 2 + ^3 + • • • (4.7) 
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We shall write down here the first three terms of the expansion, which after some tedious 
computations are given by (the expressions below hold at the quantum level as well): 



H 2 = - 




jj^n,n-l jj^n 



n 3 



+ In,n+lN„-lX, 




n-1 




where we define 




■n 



n,n—l 



= x, 



n— 1 Tn- 



(4.9) 



It is clear that as we consider higher orders in the expansion, the terms associated to the 
defect become less and less local. And although the defect is attached to a particular site 
n, its effect to higher integrals of motion becomes highly non-local. A similar behavior 
is naturally expected when deriving the relevant Lax pairs as will be transparent in the 
subsequent section. 

4.2 The associated Lax pair 

In this case one has to distinguish various cases due to the presence of the impurity. More 
precisely as we consider higher order expressions we need to take into account more and 
more points around the defect in order to include all the possible interactions. For instance, 
to derive A^ 2 ) we consider the "bulk" points and separately the point n, n + 1. For A^ 3 ) we 
separately evaluate the operator for the points n, nil, n + 2 an so on. Indeed, the main 
observation is that the presence of the defect described by the Lax operator L n induces non- 
trivial "boundary" type effects onto the neighboring operators Aj around the defect point. 
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More precisely, the generic expression for Aj the sites n, n + 1 or instance are given as: 
An (A, fi) = Ln-i(A)...Li(A) L N (X)...L n (X) 

A n+1 (A, n) = ^ L n (X)...L 1 (X) L N {X)...L n+1 (X) (4.10) 
X — fi 

and so on for points around the defect. The non-trivial "boundary" effects are due to the 
fact that the L operator is located near or on the edges of the sequence of the Lax operators 
in the latter expressions. 

After some quite tedious computations we conclude that: the Lax pair A^ 1 -* remains the 
same as in (I3.19P for all sites, A^ for j ^ n, n + 1 is given by expression (I3.19p . whereas 



H (3 n + X n+ i \ A ( 2 ) / /i Xn+l 

Also A,- for j ^ n, nil, n + 2 is given by (13.191) and 



A (2) = IH Pn + X n+l A ( 2 ) = fl X n+l 

' ' -A„_! 0i' n+1 Un-^n-l ' 



(3) / " + In-lAn-2 

-fiX n _ 2 — X n _ 3 + N n _ 2 ^n-2 — X n _2 ■En—l 



Al 3 ) 



i (3) 

Wl 



A* 2 + ^n.re+lArt-l (XXn^n+i + I n +l — N n +iX n -(-i + f 

"^A'n-l — A n _2 + N n _lA"„ A^ n _i 

/i 2 + I n +iA„ in _i flX n+ i + X n+ 2 ~~ 

,n— l^n+1 



(3) [ /i 2 + ^n+2A„ + i ^X n+ 2 + ^n+3 — N n+2 £ n + 



2 



a n+2 



'^X n+ i — X nn _i + N n+ iX n+1 — X n+ ix n+ 2 



(4.12) 



where we define 



f — Xn+2 — Xn+l ~ ®n((3 n + 2x n+ i) 

g = X n -i - X„_ 2 - «„(7n - 2X n _i). (4.13) 

Having been able to explicitly derive the first integrals of motion as well as the associated 
Lax pairs we may now identify the corresponding difference equations of motion, and check 
the consistency of the approaches followed. Indeed, both descriptions, i.e. the Hamiltonian 
as well as the zero curvature condition provide as expected the same equations of motion. 
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Let us now focus on the third charge, and extract the relevant equations of motion. These 
are given for j ^ n, nil, n ± 2 by equations (I3.2ip . whereas for the points around the 
impurity we obtain: 

N n _l + x n _ 2 N n _ 2 + X n _ 3 x n _ 2 + X n -\ 

X n -2 = — X n _4 + 2X n _3N n _2 + X„_3N n _3 — X n -2^ n -2 ~~ x n-lX 2 _ 2 ~ X n -3 

Xn-1 — x n+l ~ ^ x n,n+l^n-l — Nn+iX^+i + X n -i'N^ l _ 1 + X^_iX n —2 + x n,n+l + f 

X n -l = —X n _ 3 + 2X n _ 2 N n _i + X„_ 2 N n _ 2 — ~ \n+l^n-l — Xn-2 

x n+l — x n+3 ~ 2x n+2 N n+ i — X n+ 2N n+ 2 + X n+l^n+l + X^ +1 X ntn _i + X n+ 2 

X n+ \ = — X n ^i + 2X nj „_iN„ + i — X„_iN n _i — X n+ iN^ +1 — x n+ 2X 2 +1 — X n , n _\ + g 
2x n+3 N 

n+2 N n+3 + X n+2 K+2 + X n+2^V+1 + x n+3 

X n+ 2 = —X ntn ^i + 2X n+ iN n+ 2 + X n+ iN n+ \ — X n+ 2^ n +2 ~ X n+3X^ + 2 ~ X n+ i. 

(4.14) 

Particular attention is given to the defect point. In this case one has to take into account 
the defect degrees of freedom and the exchange relations among the elements a, f3, 7, 5 
when considering the equations of motion from the Hamiltonian. From the zero curvature 
condition on the other hand one has to bear in mind that exactly on the defect point the 
L-operator is modified to L, thus the condition may be rewritten as: 

L n (X) = A n+1 (A) L n (X) - L n (X) A n (A) (4.15) 

and the entailed equations of motion for the defect point are given as: 

«n = -/3„N n -iX n -i - 7„x n+ iN n+ i + finX n -2 + lnXn+2 - a n fi n X n _i - a n -y n x n+1 

/3 n = 2a n x n+1 N n+ i - 2a n x n+2 + 2(3 n x n+1 X n _ 1 + (3 2 n X n _i - f3 n ^ n x n+ i + 2a 2 n x n+ i + a 2 n f3 n 

7„ = 2a n X n _iN„,_i - 2a n X n ^ 2 - 2^ n x n+ iX n ^i - j n f3„X n ^ 1 + 2a 2 n X n _ x + 7^+1 - a 2 n j n . 

Taking the continuum limit of the discrete model under study is a significant aspect of 
the whole process. It is an essential step towards understanding how integrability can be 
preserved in the continuum case. There is a discussion on the continuum NLS models in [8], 
but there is no convincing argument as far as we can understand on the issue of integrability. 
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Both descriptions i.e. the Hamiltonian versus the Lax pair formulation are needed in order 
to obtain a complete view of the problem at hand. It is technically more convenient in many 
cases to use the information from the Lax pair formulation or vise versa, however in most 
cases combination of both descriptions helps to completely describe the problem especially 
when dealing with the continuum version of a lattice integrable model. 



5 The continuum limit: a first glance 

In order to proceed with the continuum limit of the discrete NLS model let us first introduce 
the spacing parameter A in the L-matrix of the discrete NLS model as well as in the L 
matrix of the defect (index free notation): 

/l + AA- A 2 xX Ax\ 
= I -AX 1 I <"> 




L(X) — AX + ( H (5.2) 



where we now define: 



a = -S = - cos(2A6>), =- sin(2A#)e 2 ^, 7 = - sin(2A6>)e- 2 ^, (5.3) 
2 2 2 



we also define: 



= Vi e -«* = Y, (5.4) 



the latter identifications will be used in the following analysis. Notice that the spectral 
parameter A is also suitably renormalized to AA in both L and L matrices in order to 
formulate a sensible continuum limit process compatible also with the continuum linear 
algebra (see also [H]). Moreover, such a renormalization is necessary if we wish the whole 
process to be compatible with the so called "power counting" argument introduced in [19|. 

Having introduced the appropriate spacing parameter in the L-operators above we may 
now consider the continuum limit of the integrals of motion of discrete NLS models and 
the associated Lax pairs. Before obtaining the continuum limit let us first introduce the 
following notation. In particular, we set: 

Xj x~(x), Xj — > X~(x), 1 < j '< n — 1, x G (—00, xo) 
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Xj -> x + (x), Xj^-X + (x) 1 n + l<j<N, x e (x , oo). (5.5) 

where Xq is the defect position in the continuum theory. Note also that in order to perform 
the continuum limit we bear in mind that: 

n-l 



A ~+ / dxf-(x) 

A E ^ ^ / dxf+(x). (5.6) 



j=n+l x o 

The continuum limit of the first integral of motion is then given as: 



H {1) = - dx x-{x)X-(x) - j dx x + (x)X + (x). (5.7) 



Notice that in the first integral we considered terms proportional to A, whereas in the 
second integral the first non trivial contribution to the continuum limit is of order A 2 . The 
respective continuum quantity reads then as 

r x o r°° 
U {2) = - / dx x-\x)X-(x) - / dx x + \x)X + (x) 



+ x (x )X (x ) - x + (x )X (x ) + x + (x )Y(xq) - y(x )X (x ) + ^y(x ) 

(5.8) 

the prime denotes derivative with respect to x. 
Note that in the continuum limij 2 ]: 

L(A)~I + AU(A) (5.10) 

also the continuum zero curvature condition with Lax pair U, V takes the from: 



U- V + 



U, V 



0. (5.11) 



2 Notice that the L matrix in the continuum limit may be expressed as: 

L - cr 3 + A U (5.9) 

(7 3 the familiar Pauli matrix. We could have chosen instead L — a^L, which also satisfies the quadratic 
relation (|3.3[) . and has the expected continuum behavior (|5.10[) . Such a choice would slightly modify the 
defect terms in the local integrals of motion. Note that such modifications can be suitably implemented in 
the continuum monodromy matrix, but we shall discuss this matter in detail elsewhere. 
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The Lax pair associated to the first integral in quite trivial and coincides with the one in 
(I3.19p . The Lax pair associated to the second integral of motion is given by the following 
expressions: 

V (2) (^, x) = 
V (2) Qu, x) = 
V (2) (^, x ) = 

Due to continuity requirements at the points also a similar argument in [18j, we 

end up with the following sewing conditions associated to the defect point: 

y(x ) = x~(x ) - x + (x ), 

Y(x ) = X-(x )-X + (x ). (5.13) 

Notice that the continuity argument may be successfully applied to the points around the 
defect, however as expected a discontinuity (jump) is observed exactly on the defect point. 
It should be emphasized that the L operator is altered at x (i.e. L — > L), leading to 
modification or discontinuity in the zero curvature condition at Xq, which accordingly lead 
to adjustments in the induced equations of motion. 

It is also quite straightforward to show that if the sewing conditions (15.131) are valid then 

{Hu U 2 } = 0, (5.14) 

which is a first good indication of the preservation of the integrability in the continuum case 
as well. However, this is somehow "on shell" information, given that one requires (15.131) in 
order to prove the Poisson commutativity (I5.14p . Moreover, the constraints (I5.13P provide a 
first hint on the existence of an underlying non-ultra local algebra (see e.g. [20J) associated 
to the defect point. And although in the discrete case one deals with an ultra local algebra 
there is an indication that in the continuum limit one has to consider a generalized non-ultra 
local algebra in order to efficiently describe the point like defect at xq. 

Here, we only provide a first glimpse on the continuum limit of the discrete NLS models. 
To obtain the continuum counterparts of the higher integrals of motion and the associated 



/i x (x) 
-X~(x) 

(I x + (x) 
-X + (x) 

fi x + (x ) + y(x ) 
-X-(xo) 

/i x + (x ) 
Y(x )-X-(x ) 



x G (— oo, x Q 
x E (xq , oo) 



(5.12) 
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Lax pairs requires subtle manipulations. Such an explicit construction is beyond the intended 
scope of the present investigation, however we shall analyze this intriguing issue in full detail 
in future works [2T] . 

6 Discussion 

Let us briefly summarize the main findings of the present study. The main aim of this work 
was the investigation of the discrete NLS model in the present of an ultra local integrable 
defect. Based on purely algebraic considerations we were able to extract the first charges 
in involution for the discrete model. The model is by construction integrable given that the 
bulk L matrices as well as the Lax matrix associated to the defect are required to obey the 
same Poisson bracket structure. 

Then, again by exploiting the underlying algebra, we extracted the associated Lax pairs. 
Particular attention was given to the construction of the Lax pairs around the defect point. 
It turned out that the behavior of more and more points around the defect is affected as 
we move to higher order expressions. Having this information at our disposal we were able 
to derive the sets of the difference equations of motion. Finally, we provided a first insight 
of the continuum behavior of the system by considering the continuum limit of the first 
two integrals of motion. This led to certain sewing or compatibility conditions that ensure 
Poisson commutativity of the first two integrals of motion at the continuum well. 

It is worth noting that at this stage it is difficult to conclude whether or not our results 
correspond to the results of e.g. [S] or if any comparison whatsoever can be made, given 
that the issue of integrability is still open in [5] . The systematic continuum limit of discrete 
models in the presence of integrable defects in the spirit of [19] is our next target |21j . and 
it will provide a deeper understanding on the connection with earlier works. In fact, this 
systematic process will lead to the derivation of the continuum limit of higher integrals of 
motion, such as T-L^ and the corresponding Lax pairs, as well as the associated constraints 
(sewing conditions). Note that the entailed higher constraints will involve as expected spa- 
tial derivatives. Moreover, Poisson commutativity of all the entailed charges needs to be 
explicitly checked so that we can claim that integrability holds. Compatibility of the higher 
sewing conditions should be also explicitly checked. These are highly non-trivial technical 
points, and will be presented in full detail in future investigations. 
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Similar ideas may be put forward in the case of other well known prototype models such 
as the Heisenberg model (see e.g. |TT]) aiming also at the investigation of the corresponding 
continuum theories. In any case, a detailed analysis on continuum integrable models in 
the presence of defects turns out to be a fundamental issue, which will be addressed in 
forthcoming publications. 
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